Dispersive Charge and Flux Qubit Readout as a 
Quantum Measurement Process 

Lars Tornberg and Goran Johansson 

Applied Quantum Physics, MC2, Chalmers, 96 Goteborg, Sweden 



We analyze the dispersive readout of superconducting charge and flux qubits 
as a quantum measurement process. The measurement oscillator frequency 
is considered much lower than the qubit frequency. This regime is interesting 
because large detuning allows for strong coupling between the measurement 
oscillator and the signal transmission line, thus allowing for fast readout. 
Due to the large detuning we may not use the rotating wave approximation in 
the oscillator- qubit coupling. Instead we start from an approximation where 
the qubit follows the oscillator adiabatically, and show that non-adiabatic 
corrections are small. We find analytic expressions for the measurement 
time, as well as for the back-action, both while measuring and in the off- 
state. The quantum efficiency is found to be unity within our approximation, 
both for charge and flux qubit readout. 



PACS numbers: 74.50.+r 03.65.-w 03.67.Lx 42.50.Lc 



1. INTRODUCTION 

Superconducting qubits are considered one of the technologies capable 
of bringing us quantum information processing (QIP) devices in the future. 1 
One cornerstone of QIP is the possibility to convert the quantum state into 
classical information, i.e. read-out. A large effort has gone into optimizing 
different read-out techniques, see e.g. Ref. 2 and references therein. Of 
course reading out a qubit only makes sense if it can stay coherent long 
enough to contain true quantum information. Being solid-state qubits, su- 
perconducting qubits interact rather strongly with their environment, mak- 
ing dephasing the worst enemy of superconducting QIP. 
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Quite naturally, which type of noise is most dangerous depends on the 
qubit design. For example charge qubits 3 are most sensitive to charge fluctu- 
ations, 4 while flux qubits are most sensitive to flux noise 5 and charge-phase 
qubits 6 are sensitive to both. 7 At the degeneracy point the two charge/flux 
qubit states have equal average charge/flux. Thus, at this point they are 
protected against environmental charge/flux noise to first order in the cou- 
pling. To minimize the effect of the typical noise, qubits are now operated 
so that they remain at this degeneracy point as much as possible. 

The natural readout for a charge/flux qubit is to measure the aver- 
age charge/flux of the qubit. Since the qubit states have equal average 
charge/flux at the degeneracy point, the natural type of readout is no longer 
possible. This initiated an interest in dispersive readout techniques, where 
the qubit state influences the properties of a harmonic oscillator, which is 
then probed. 8-10 

In the experiment of Wallraff et. al. at Yale, 10 a charge qubit was 
coupled to an oscillator, realized by a coplanar waveguide. Using this setup 
they managed to determine the state of the qubit, while maintaining long 
coherence times. In this case, the oscillator and qubit where nearly degen- 
erate in frequency and an oscillator with a high quality factor Q was used 
to shield the qubit from relaxation. 

In this paper, we concentrate on the situation where the qubit and os- 
cillator frequencies differ by orders of magnitude. This will allow us to use a 
low Q oscillator, thereby increasing the read-out speed of the circuit. Here it 
is the detuning that shields the qubit from measurement induced mixing, and 
also allows relaxation to be minimized through effective filtering of the trans- 
mission line at the qubit frequency. In this setup the state of a charge/phase 
qubit can be read out through its effective capacitance/inductance, which 
shifts the resonant frequency. At the degeneracy point, the state dependent 
quantum capacitance/inductance for the two states have the same magni- 
tude but differ in sign. For the charge qubit the quantum capacitance was 
recently measured by two groups. 11 ' 12 

The theory of quantum measurement states the impossibility to discern 
the state of a quantum system without destroying the phase coherence in 
the system. This promotes the idea of defining the quantum efficiency rj of a 
measurement as the ratio between the dephasing time tu, and measurement 
time t rns , bounded above by rj = t ip /2t ms < I. 13 (To be consistent with 
the results of the Yale group, 16 we have decided to use their definition of 
measurement time. This, together with the desire to keep rj < 1 is the reason 
for the appearance of the one half in the definition) . The quantum efficiency 
gives a measure on the back-action of the read-out and deriving 77 in terms of 
circuit parameters is thus a good starting point for optimizing the read-out. 
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In Ref. 14, together with C. M. Wilson, we investigated the quantum 
efficiency of the read-out proposal mentioned above, as well as discussed the 
performance of a realistic circuit design. Here we give the more detailed the- 
oretical background for those results. The large detuning between the qubit 
and oscillator makes standard methods of quantum optics 15 ' 16 inappropri- 
ate. Instead, we start from an approximation where the qubit follows the 
oscillator adiabatically, and show that non-adiabatic corrections are small. 
Our results are based on the quantum network theory, introduced by Yurke 
and Denker. 17 As already discussed in Ref. 14, the efficiency for both charge 
and phase qubit readout is ideal within our approximation. This result is 
independent of the oscillator Q-value. Furthermore, we find that a low Q 
oscillator shields the qubit from thermal dephasing. Finally we treat explic- 
itly the limiting case of a transmission line capacitively coupled directly to 
a charge qubit. Being of more academic interest, we find an ideal quantum 
efficiency also in this case. 



2. CHARGE QUBIT 

The goal of this section is to derive the Hamiltonian for the super- 
conducting charge qubit complete with read-out device. Having the full 
quantum description of the circuit we then make suitable approximations 
to derive the quantum capacitance of the Cooper-pair box. This gives us 
an adiabatic Hamiltonian which we can solve to obtain the dynamics of the 
system exactly. In sectional we give bounds for the non-adiabatic contri- 
butions, and justify our approximations in detail. In the rest of the paper 
we use the adiabatic Hamiltonian to investigate the efficiency of our mea- 
surement scheme. 

The single Cooper-pair box operated as a charge qubit 18 can be viewed in 
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Fig. 1. Cooper-pair box operated as a charge-qubit complete with read-out 
circuitry, see text. 



L. Tornberg and G. Johansson 



Fig. ^ A small superonducting island is connected to a superconducting 
lead through a Josephson junction with Josephson energy Ej and capac- 
itance Cj. The island is also coupled to a gate voltage V g via the gate 
capacitance C g . By tuning V g , the desired working point of the qubit can 
be reached. To enable readout an LC-oscillator is coupled to the box via a 
measuring capacitance C m . The oscillator is in turn coupled to a transmis- 
sion line, here modelled as a semi-infinite line of LC-oscillator s in parallel. 
The line is characterized by its capacitance Ct and inductance Lt per unit 
length. Through this line, all measurements of the qubit will be performed. 

We choose, as our starting point, to write down the Lagrangian for the 
system. 17 ' 19 This method gives automatically the conjugated variables of the 
system and thus provides us with the dynamics. Since the circuit contains 
a Josephson element the natural choice of coordinates is the phase across 
each circuit element, which is related to the voltage drop by 4>i = I Vidt. 
Following standard procedure 19 the capacitive and inductive energy of the 
circuit are the kinetic and potential energy of the Lagrangian 

2 2 2 2 2 ^2 

where <i>o = h/2e is the flux quantum. Applying Kirchoff's voltage law 
gives the constraints on the circuit coordinates which gives the Lagrangian, 
L = T — V , for the system 

C qb &j (Cqsc + CJ&c C c {^f $| /gj 
2 2 2 2L J V$o 

- C m $ c *J ~ C c ^ P i ~ CgVgZj + + 

where we have introduced the capacitances C q b = Cj + C g + C m , and 
C sc = C + C m for brevity. To simplify the derivation of the Hamiltonian 
we introduce the vector notation for our coordinates 



$j $ c §1 ... . (3) 



The Lagrangian can now be written more compactly 
L _ -* « - C a Vfij - -f L + Ej c„ 8 ) - £ A,Lg^, (4) 
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where C is the mass matrix of the system. This can be divided into C sys , 
describing the capacitative elements in the qubit-oscillator circuit and a pure 
diagonal matrix Ctl which only contains information about the capacitance 
per unit length in the transmission line 



c = 


Csys 








Ctl 




Ct 


.. 


Ctl = 





Ct 
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sys 



Cqb —C m 

Cm C sc ~t~ C c C c 

-C c C c + C T Ax 



(5) 



The Lagrangian contains all information needed for deriving the system dy- 
namics. However, when doing quantum mechanical calculations one often 
prefers to work with the Hamiltonian. We thus Legendre transform 20 the 



Lagrangian following Appendix A. and arrive at 

H = H q b + H osc + Hi n t + Htl, 



(6) 



where H q b describes the Josephson junction alone and also the interaction 
between JJ and oscillator 



H gb = ^(PJ + CgVg f + iDnpc + D nPp ) (pj+C g V g )-E j cos (^f 



(7) 



and where Dij are the elements of the matrix 



n = n- 1 = 1 

^sys n ,r — 
^qb^osc 



Cose C m 
Cm Cqb 
C m Cqb (CqbC 



C m 
Cqb 

C m ) I c c 



a 



qb 



(8) 



The last three terms in the Hamiltonian, Eq. @, describe the oscillator, the 
transmission line and the coupling between them; 



Htl 



d 22 , e 

— Pc+ 2L 



D 



1 



33 2 , 




(9) 



where the operators pc,PJ,Pp and p\ are the conjugate momenta to $c", $ j, <3?f 
and <5f respectively, and have dimension of charge. By definition, the 
phases and charges (pi and pj obey the canonical commutations relation 
= ifidij. The kinetic energy of the Hamiltonian is thus represented 
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by the the charging energies of the circuit. The D-matrix contains informa- 
tion about the effective capacitances associated with the respective charges. 

The J J island is caracterized by a charging energy, Ec = Due 2 / 2 as 
e 2 /2C q b and the Josephson energy Ej. To realize a two level system, the 
charging energy of the island needs to be much larger than the Josephson 
energy. For this choice of parameters we can limit the number of excess 
electrons on the island to the values {0, 2}. Projecting the qubit Hamiltonian 
on this subspace of the whole Hilbert space in charge basis, we get the usual 
qubit Hamiltonian in the language of Pauli spin matrices 



where k = Dx^/Dn = C m /C osc is a dimensionless coupling constant and 
E e i = 4E C (1 — n) is the electrostatic energy of the island depending on both 
gate charge and induced charge from the oscillator and transmission line 
n = uq + k(tlc + n p ). The number operators n are related to the charges by 
n c = Pc/e, n p = Pp/e and no = C g V g /e respectively. The last term in Eq. 
(flfl|) is just an offset in the energy of the system and thus discarded. 



During read-out the oscillator will be driven close to the bare resonance 
frequency which we consider much smaller than the qubit level-splitting 
Ej/h. This difference in frequency allows us to make an adiabatic approxi- 
mation, separating the qubit and oscillator dynamics. Thus the qubit eigen- 
values and eigenstates will depend on the oscillator and transmission line 
coordinates. However the levels will never cross and the qubit is frozen in 
its state. In this regime, the qubit can be described by an effective adiabatic 
Hamiltonian 



which is obtained by a unitary transformation of the the qubit Hamiltonian 
in Eq. (jlOft . For more details on this, see Sec. E]The expression in Eq. 
is similar to the usual expression for the Hamiltonian of a charge qubit in 
computational basis. The difference in this case is that the working point of 
the qubit is determined not only by the V g dc-bias, but also on the induced 
charge on the oscillator. The second term is usually proportional to the 
identity operator and ignored. Here it contains both system and transmission 
line operators, and can not be thrown away. However, it only represents a 




(10) 



2.1. QUANTUM CAPACITANCE 




(11) 
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small shift in pc and p\ and can thus be included into a new set of shifted 
operators. Considering a small coupling between qubit and oscillator we can 
expand the qubit energy to second order in 5n = n{nc + n p ). The effect 
of the qubit can thus be summarized in three terms; one state dependent 
off-set which does not affect the circuit dynamics, a linear term in k{uc + 
n p ) which is zero at the degeneracy point, and finally a term quadratic in 
oscillator and transmission line charge. Discarding the two first we arrive at 
the Hamiltonian determining the low frequency dynamics of the oscillator 
and transmission line, at the charge degeneracy point 

H = ^ + § + f^c+P P ) 2 

+ — p P + S ? {-^C7~ + 2Tt J + D23PCPp - (12) 

Here we have defined the coupling between the charge qubit and oscillator 
gc = Cq/Cq SC = 8 *z c e , and the quantum capacitance as 

^-e^T- (13) 

The motivation behind this definition comes from the fact that the qubit 
shifts the electrostatic energy of the oscillator, which can be interpreted as a 
shift of the oscillator capacitance, C osc — ► C osc + Cq. This interpretation is 
possible as long as Cq <C C osc . The corresponding qubit induced shift of the 

oscillators resonance frequency is simply given by u osc = 1/ J L(C osc + Cq). 



3. FLUX QUBIT AND QUANTUM INDUCTANCE 

The radio-frequency superconducting quantum interference device (rf- 
SQUID) operated as a flux qubit 21 can be seen in Fig. [21 The readout part of 
our analysis also applies to the 3-junction persistent current qubit, 5 but the 
rf-SQUID is chosen for simplicity. In the previous section the charge states of 
the single Cooper-pair box were used as the computational basis of the qubit. 
Here the qubit loop is biased by an external flux $> x , and the |0) and |1) 
states are realized by the clock and anticlockwise circulating supercurrents 
in the loop. The read-out circuit is the same as for the charge qubit. It is 
however inductively coupled to the qubit via the mutual inductance M, as 
opposed to the capacitative coupling in the charge case. The derivation of 
the Hamiltonian for this circuit is almost analogous to the analysis for the 
charge qubit. Making the adiabatic approximation and rotating to the qubit 
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Fig. 2. Radio-frequency SQUID operated as a flux-qubit, complete with 
read-out circuitry, see text. 

eigenbasis yields the following adiabatic Hamiltonian for the oscillator and 
transmission line, considering the qubit biased at its degeneracy point 

' , , • M2 (<t>) 2 A**, (PpT , PpP , 




($ p , - $ p ) 2 \ 

+ ^i^gf+^V" 1 )- (14) 

where C c = (C + Ci n )/(CCi n ), (cp) is the magnitude of the flux generated 
by the qubit persistent current in a current eigenstate, and A is the qubit 
energy splitting. Just as in the charge case, the slow transverse coupling to 
the oscillator is seen as a longitudinal coupling in the eigenbasis of the qubit. 
Due to the phase-coupling between oscillator and qubit we get the coupling 
constant gi = Lq/L 2 where the quantum inductance is defined as 

Lq ~ (15) 

so that the coupling can be incorporated into an effective inductance L g / e = 
L± Lq. 



4. NON- ADIABATIC CORRECTIONS AND QUBIT 

RELAXATION 

In the following we will calculate the dynamics of the oscillator and 
transmission line, considering the qubit state to be fixed in its adiabatic 
eigenbasis. This is of course an approximation, and there are two different 
effects induced by the measurement circuit that will act to change the qubit 
state. First, too strong driving will induce transitions, which in the large 
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amplitude limit will be of Landau-Zener type. Secondly, seen from the qubit 
perspective, the measurement circuit consists of a dissipative transmission 
line filtered through the oscillator. This will add to the relaxation rate of 
the qubit, also when not measuring. 



4.1. Non-Adiabatic Corrections 



Equation Ijllj) was derived by diagonalizing the qubit Hamiltonian, a 
transformation which depends on the state of both oscillator and transmis- 
sion line. In the adiabatic regime however, it is possible to find eigenvalues 
and eigenstates of the qubit for each value of the charge on the oscillator 
and transmission line. Because of the above, the transformation will how- 
ever also affect the oscillator part of the Hamiltonian, and we must consider 
this effect and make sure that it is negligible. We begin with the unrotated 
Hamiltonian for qubit and oscillator (see Eq. (fTT 



H 



E e i(n) 



Ej , p 

o T A 

2 2C 



2 $2 
+ 



2L 



(16) 



where the coupling is in E e i, as assumed before. Making the unitary trans- 
formation to qubit eigenbasis 22 



H 



fHu, 



cos rj/2 
sin 77/2 



— sin 77/2 
cos 77/2 



77 = arccot 



E e i{n) 
Ej 



(17) 

diagonalizes the qubit part of H, but also affects the oscillator term that 
contains the phase. At the degeneracy point no = 1 the transformation is 
given to first order in k 



u ~ 



V2 




(18) 



where x = 4\E c Kn/Ej. The term that we have disregarded when deriving 
Eq. (JTTJ) is thus 



.$2 

u { — u 
2L 




ih 



1 
-1 



(19) 



which would induce transitions in the qubit. It must therefore be compared 
with the natural transition energy in the qubit, given by the level splitting 
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Ej. Since we explicitly assume both small amplitude Ak5uE c <C Ej and low 
frequency %uj <C Ej of the driving, the ratio is given by 

4_nE^SnJ_ <<; ^ (2Q) 
-Ej Ej 

which is indeed negligible within our approximation. We note the similarity 
of the above condition and the one needed for neglecting Landau-Zener(LZ) 
transitions. 23 The LZ transition rate is given by 

T LZ oc e - 2 ^y n \ v = d t E eh (21) 

which thus also can be neglected if the condition in Eq. Q2U[) is fulfilled. For 
the flux qubit the corresponding condition for small amplitude driving is 

M\(6)Vk , , 

~kljt « L (22) 

In Sec. ISHwe will find that the quantum efficiency is independent of the drive 
strength. But the speed of the measurement is proportional to the power 
of the drive, as we will see in Sec ED Thus, in a real experiment one would 
like to drive as strongly as possible, but one should then make sure not to 
violate Eos. (1271 or (|22l . 



4.2. Qubit Relaxation 

Although the pure dephasing induced by the transversely coupled oscil- 
lator circuit vanishes to first order in the oscillator-qubit coupling, the qubit 
relaxation rate Ti has a first order contribution. Thus we evaluate the relax- 
ation rate using standard weak coupling expressions, 22 giving that it is pro- 
portional to the real part of the impedance seen from the qubit, at the qubit 
frequency uj q f,. For the simple oscillator circuits shown in Figs.^ an d EJwith 
a low Q, the relaxation induced by the transmission line can easily dominate 
over other relaxation mechanisms. The solution to this apparent problem 
is to add a non-dissipative low-pass filter between the transmission line and 
the oscillator. Due to the large frequency difference between the oscillator 
and qubit, this is straightforward. The high-frequency electromagnetic en- 
vironment responsible for qubit relaxation can be designed independently of 
the low frequency environment which should let the readout signal through. 
In principle, a commercially available pi-filter can improve the relaxation 
rate by a factor of 1000, while not affecting the low frequency measurement 
properties. For realistic designs 14 this is well beyond the point where other 
sources of relaxation will dominate. 
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5. SYSTEM DYNAMICS 

The questions that we address in this paper concern the efficiency of 
read-out by probing the quantum capacitance (inductance) of a charge (flux) 
qubit. The basic idea is that we send a pulse through the transmission line, 
which will be reflected by the combined oscillator-qubit circuit. In our ap- 
proximation the pulse will see two different harmonic oscillators, depending 
on the qubit state. While the information about the state of the qubit state 
thus is transferred to the reflected pulse, the unavoidable photon number 
fluctuations in the pulse will dephase the qubit. The statistics of these fluc- 
tuations is determined by the way the pulse is generated, and the dynamics 
of the transmission line coupled to a harmonic oscillator. We thus start out 
by calculating the dynamics of the free transmission line. Given this, the 
qubit-oscillator circuit will act as a boundary condition to this solution. 



5.1. Equations of Motion: Transmission Line 



The time-evolution for the transmission line degrees of freedom is gov- 
erned by the Htl in Eq. ©• For x > we consider the transmission 
line phases and charges in the limit Ax — > 0. In this limit the Heisenberg 
equations yield the coupled equations, two for each i > 



Pi 



Pi 



Ax (gj+i - $j) 2 - (gj - ^_i) 2 
(Ax) 2 



dx d 2 ®i 
Lt dx 2 ' 



2L T 

Pl = C T dx<&i. (23) 

Combining these two equations to one we see that the phase in the trans- 
mission lines obeys the massless scalar Klein-Gordon equation 

<9 2 $ 1 <9 2 $ 



dt 2 



(24) 



LtCt dx 2 

which indicate the bosonic nature of the transmission line excitations. Equa- 
tion (|24[) has a solution in the form of right and leftgoing travelling waves 
which formally can be written as 



$ = $ m - + t + $ 



\out 



X 

— +t 

V 



(25) 



where v = 1/ 1 \JLtCt is the velocity of the waves in the transmission line. 
By differentiating Eq. (|25|) we can get a relation between the partial deriva- 
tives 17 



1 <93> 
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The spatial derivative of the probing field at x = together with Eq. 
will later provide us with the coupling between circuit and transmission line. 
It is however natural to first derive the statistics of the transmission line and 
then go on to the circuit dynamics. 



5.2. Expressing the Transmission Line Using Creation and 

Annihilation Operators 

The Klein-Gordon equation of motion is easily solved by introducing 
the Fourier transformed operators in the transmission line 

1 POO 

1 I ~ » • - * ikx . 



*(s,t) = -=/ <s>(k,ty« x dk, (27) 

V27T J-00 

and we get an ordinary differential equation for each value of the wave- 
number k. Solving Eq. (|24jl is reduced to the problem of solving 

$(fc,t) + v 2 |fc| 2 *(M) = 0, (28) 

which is the classical equation of motion for a harmonic oscillator with the 
dispersion relation uJk = v\k\. In analogy with the ordinary harmonic oscil- 
lator it is convenient to work with the creation and annihilation operators 
a and a) instead of charges and phases. However, since we are dealing with 
a field equation each Fourier component is treated independently and we 
write 24 



*(M) = ^2^7 k ( ak(t) + a - k{t) )> (29) 

where a and a* obey the canonical commutation relations [a&, a',,] = 5(k — k') 
and [afc, ay] = 0, which give the statistics of the transmission line excitations 
(photons). The harmonic oscillator nature of the transmission line makes 
the phase of the creation and annihilation operators rotate with angular 
frequency u^, a^(t) = aj t e~ luJkt , at(i) = a\,e luJkt and we arrive at the final 
expression for $ in terms of a& and at 



/ 4ttC t J 


f°° dk , 




-oo \A^k 



$(x,t) = J-r^r I (afce-*^**-^ + 4e iK *" fex) ) • (30) 



In our calculations, we use the in-(out) formalism introduced in Eq. (|25[). 
to derive the circuit scattering matrix. For this purpose it is convenient to 
split the integral in Eq. ()30|) into positive and negative wavevectors, which 
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at x = can be recognized as the in and out-fields of Eq. (|25|). As a final 
remark, one may also write the fields in terms of integrals over frequency 
giving the normalization 



$(x = 0,t) 



47T 



a u e 



-iuit 



(31) 



where the annihilation and creation operators in frequency space obey the 
canonical commutation relations [a^,a^,] = 5(u> — a/) and [c^cv] = 0. 



5.3. Equations of Motion: Charge Qubit 



Considering the adiabatic Hamiltoninan in Eq. (|12|) and the relation in 
Eq. (|26j) . the Heisenberg equations of motion for the circuit are 

= (£>23 + 9°z){pc +P P ), ®p = (D 33 + ga z )p p + (D 2S + gcr z )p. 

(32) 

This set of four coupled equations can easily be reduced to two, only contain- 
ing the phase operators of the oscillator and the transmission line. Since the 
effective circuit only contains linear elements it is straightforward to solve 
Eq. 1)32(1 in Fourier space. The linearity is a direct consequence of the weak 
coupling between the oscillator and the qubit, where we have replaced the 
non-linear Josephson element with a linear capacitance. Both the outgoing 
field and the charge operator coupling to the qubit are expressed in terms 
of the incoming field 

. 9? T i 

= -^$ m ( W ) = X9/e H$ m ( W ), 

*&» = ^^^) = e^^\ (33) 



where Zq = Lt JCt is the characteristic impedance of the transmission 
line. We also introduce the short-hand notation D^ e = Di + ga z . Since 
there is no dissipation in the lumped circuit we have |$ (u)\ = |$ m (u;)| 
and the in- and outgoing field only differ by a state-dependent phase. In the 
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language of Subs. 15.2.1 we write Eqs. (|33|) in the time domain 



47T Jo 



>hZ n 



du> 



^ I ^( X9/e ( W )fl m Me-^)+h.c.), (34) 



where the operators contributing to the charge on the oscillator are written 
as one operator p g / e (t) = p^ e (t) + Pp^ e (t). We now have the exact dynamics 
for the combined circuit + transmission line system. However, to get simple 
analytical expressions for the measurement and dephasing time we need to 
write the relations N(uj) and x( w ) on a more simple form. We thus make a 
Breit-Wigner approximation, see | Appendix B. Thus we write as 



<?/ e 



1 + i2Q9l e 



,s/ e 



,s/ e 



g/ e 



23 



\L{l + C c Dli e ) 



Q 



9/e 



ZoClL^tY 



(35) 



where the peak value Xo ^ s evaluated at the resonance frequency uj 
and Q g l e is the quality factor of the oscillator. 



' 



5.4. Equations of Motion: Flux Qubit 

The flux-qubit is treated in exact analogy with the charge qubit. The 
only major difference is that the qubit degree of freedom now couples to the 
phase, and not charge operator of the oscillator. Given the Hamiltonian in 
Eq. (|T4*|) we solve for $ and <fr° ut in terms of the in-field & l p n 

N 9 ' e {uj) = C(l-iC c Z uj)(l-L 9 / e Cuj 2 )-C c 

= w/^ ffi") = e^'^^C"). ( 36 ) 

where x( w ) f° r the flux qubit can be written on the exact same Breit-Wigner 
form as in Eq. Q35JI . The only difference lies in the parameters describing 
the peak height, width and resonance frequency. In this case they are 



/e _ 2C C Q"° 9/e _ 1 ni/e-& »!!l f37) 

Xo ~ 1 c ' "° "v^cw' Q -z c^c cS > (37) 
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where C eS = C/(l - C c /C). 



5.5. Phase of the Reflected Signal 



We now have a full quantum description of oscillator and qubit dy- 
namics in terms of the incoming field. The interaction can be written on 
a simple form and we can use this to derive simple analytical expressions 
for the quantum efficiency for the different qubits. Using the Breit-Wigner 
approximation, the reflected phase has the following simple form 



The Breit-Wigner approximation rely on the presence of a relative pro- 
nounced resonance in the function x( w )> implying an oscillator with a finite 
quality factor. We emphasize however that this need not to be the case. 
Numerical calculations, using the exact form for N(lv) and x( w ) also & ve a 
quantum efficiency r) ~ 1. In Sec. ED we calculate rj without the presence of 
the oscillator, thereby taking away the resonance altogether. Even in this 
case, we find a quantum efficiency of unity. 

6. DETECTING THE QUBIT STATE BY HOMODYNE 

DETECTION 

The interaction between qubit and oscillator allows for an indirect mea- 
surement of the qubit state. Since the interaction is purely reactive and 
no dissipation is present in the process, all information about the qubit is 
contained in the phase of the reflected signal. To analyze the quantum effi- 
ciency of such a read-out we use standard techniques from quantum optics 
and describe it as a homodyne measurement. 15 In this setup, the reflected 
signal from the circuit is mixed with a phase-locked local oscillator field us- 
ing a linear mixer, a beam-splitter combined with a photon detector in the 
optical case. By detecting the power of this mixed signal, information about 
the phase can be extracted. Labelling the mixed signal with the operator b 
gives the field intensity at the detector 




(38) 




(39) 



Detecting intensity is equivalent to detecting photons, a destructive process 
where photons are converted into current by a photoamplifier. Our signal 
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will thus be given by the annihilation operator at the detector 

b(t) = ra(t) + rv(t) + ta(t), (40) 

which contains the average a and residual v of the LO mode. The reflected 
signal from the circuit is given by the operator a. The reflection and trans- 
mission coefficients r and t are assumed to be real. For optimal efficiency 
we need that r <C 1 which allows almost all the signal from the circuit to 
be transmitted to the detector. By amplifying the local oscillator field ra{t) 
can still be made dominant so that second order contributions in v(t) and 
a(t) can be neglected in Eq. (|39j) . Both the LO-field and signal are modelled 
as a Glauber state 25 

|{aH» = exp (j <L>[a(u)(aW - <*»(#]) |0>, (41) 

where a(uj) is the Fourier-transform of our drive signal, and |0) is the con- 
tinuum vacuum field a u \0) = 0. The drive source is nearly monochromatic 
and has a narrow bandwidth distribution (r^ <C <-Od) which is Gaussian in 
frequency 

a(u)=a ^ = , (42) 



Id 

where «o is a dimensionless constant and cod is the drive frequency. At 
the detector we observe the quantum mechanical expectation value of the 
operator N. Assuming that the deviation v from the average LO signal is 
in a vacuum state, this is given by 

(N(T)) = T(r 2 a 2 LO + 2tra LO cos(0 + 9)), (43) 

where (j> and 6 are the phases of the signal and LO-field respectively. The first 
term only contains information about the LO-field and does not contribute 
to the phase information. We thus take the second term as the signal at the 
detector S. The standard deviation of the photon number operator is 

AN = ra LO VT, (44) 

which we take to be the signal noise No. The measurement time t ms is 
defined as the time for which the signal to noise ratio is one. Choosing 
9 + (<f)9 + <j%)/2 = vr/2, and t ~ 1 it is given by 

11 

tm.s — ! F a TT •> (45 J 



4 s i n 2 



<Pr-<Pr 



where Ti n is rate of photons sent down to the circuit by the drive. We now 
have a formula for the measurement time, determined by the driving strength 
and the difference in the reflected phase between the two qubit states. This 
we may compare with the measurement induced dephasing. 
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7. QUBIT DEPHASING 



The dephasing rate r„ is defined as the rate with which the qubit looses 
its phase coherence. That is how fast the off-diagonal matrix elements in the 
reduced qubit density matrix go to zero. Since the phase coherence is the 
heart of all quantum computation, the qubit is useless after a time given by 
the dephasing time The short dephasing time is the major bottleneck 

preventing the progress in solid state quantum computing. The reason for 
such short life-times are the ever present strong couplings between qubit 
and environment inducing all kinds of fluctuations in the qubit. We must 
therefore see to that all additional sources of decoherence are minimized to 
reduce this effect. One such source is the transmission line that we have 
attached to the qubit for read-out. Photons are the elementary excitations 
of energy in the transmission line and fluctuations in these will dephase the 
qubit. In this section we calculate the dephasing rate for two different cases. 
In the first case the state of the in-field is set to a Glauber state as was 
done in Sec. ED to simulate the field generated by a monochromatic source 
used for readout. In the second thermal distribution of photons is 

used to simulate the dephasing due to the small but finite temperature in 
the experimental setup, which occurs also when not reading out. 



We begin by considering a general pure state characterizing the qubit 
and oscillator 



where a and b are complex numbers fulfilling |a| 2 + \b\ 2 = 1. At the degen- 
eracy point the Hamiltonian describing the qubit and interaction with the 
oscillator is given by 



where we incorporate the average value of p into the Hamiltonian Hq. Our 
first task is to calculate the time evolution of the <r+ operator. We use the 
charge qubit as example, but both the formalism and final result apply to 
the flux qubit as well. For our purposes it is convenient to make a unitary 
transformation to the rotating frame of the unperturbed qubit system Hq 



7.1. Dephasing from Field Correlators 




(46) 




(47) 



I (t)=e-rJ H ° dt O H (t)e*f Hodt 



(48) 
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where Hq is the unperturbed Hamiltonian. In this frame the equation of 
motion for the <r+-operator is 



~dt 



h 2 Jo 



t r 



a + {t)p\t')p\t)+p\t)a + (t)p 2 {t')+ 



dt' 



(49) 



to second order in g. Moreover, the evolution of cr+ is assumed to be slow on 
the time-scale of the oscillator. The dephasing is now given by multiplying 
Eq. (|49|) by the state in Eq. (|46|) from the left and right and tracing out 
the oscillator degrees of freedom, giving an equation of motion for ab* , the 
off-diagonal element in the reduced density matrix 



d(ab*) 



-ab* 



Ah 2 Jo 



(pWJt)) + (p 2 e (t)p 2 M')) + 



+( p j(t')p 2 Jt)) + (pi(t)pi(t>)) 



dt'. 



(50) 



The index on the p 2 -operator indicates the state of the field. We now have 
an expression for the dephasing rate. The only remaining task is to put the 
in- field in a given state and calculate the two-time correlators of Eq. (|5fl|) . 



7.2. Dephasing Induced by the Read-Out 

In the read-out process the oscillator state in Eq. (|4*6|) is a coherent 
Glauber state, see Eq. (|4*T1) 

|tf> = a\g) \a{uj)) g + b\e) ® \a(u)) e . (51) 

The oscillator state is now an eigenstate of the annihilation operator which 
makes the calculation of the correlators in Eq. 1)50(1 a matter of putting 
the creation and annihilation operators in normal order. Using the narrow 
bandwidth of the drive as given in Eq. 1)42(1 the correlators are found to be 

Air Jo u> 

16ir z Jo Jo u}\ 0J2 

The dephasing rate is given by taking the zero frequency Fourier transform. 
We make the variable transformation t — t' = T and extend the upper limit of 
the integral to infinity, an approximation which is valid if we limit ourselves 
to times longer than the field correlation time. Neglecting the imaginary 
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part of the transform, which only gives a small renormalization to the qubit 
energy splitting, we get the simple expression 

* { r^WmA = r ,„^ ^!bgM£ . (53) 

where Tj n is the rate of photons sent down the transmission line and ojd is the 
drive frequency. Solving Eq. (|5()[) together with the result from above, the 
off-diagonal element will exhibit an exponential decay with the dephasing 
rate given by 



rv = r m |-^(| X9 (^)| 2 + |x e (^)| 2 ) . (54) 

This result is valid for both charge and flux qubits, choosing the correct 
expressions for coupling g = gc = Cq/C 2 sc and g = gi = Lq/L 2 , as well as 
for the transfer function Xg/ei 00 ) from Eqs. p3|) and (|36|) . As expected the 
dephasing rate increases with the drive power. We also note that the rate 
is inversely proportional to the drive frequency squared, not surprising since 
the dephasing is given by the zero-frequency fluctuations in the field. 



7.3. Finite Temperature Dephasing 



Because of the non-zero temperature in the environment, photons in the 
transmission line will be thermally excited and contribute to the dephasing 
of the qubit. It is desireble to reduce this effect by an efficient design. We are 
therefore interested to see how this dephasing scales with the different circuit 
parameters. The state of the oscillator in Eq. 1)46)1 is given by a thermal 
distribution of photons. The correlators in Eq. (|50[) are characterized by 

(ala^) = 5(u - u')n(u), n(u) = ^ _ - , (55) 

where n(uj) is the Bose-occupation number. Given this it is straightforward 
to calculate the dephasing rate as was done in Subs. 17.2.1 Following Eq. 
dSOj) we calculate the correlator 



<p£(*)j£(0> = <p£><p?> 



'I6vr 2 



duduj' x 



x«MxM") (n((j) + 1)e -Mt-f) + ^M^H) 
^yW) + l)e-^ + tiW*>W n V)e*>'W 



(56) 
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where the first term vanishes since the p 2 -operator has zero average. The 
dephasing rate is given by the Fourier transform of Eq. (|56|) . As in Subs. 
17.2.1 we make the substitution r = t — t 1 and ignore the imaginary part of 
the integral, giving 

T * = itL ( |X9M|2 + ,XeM|2 ) (nH + 1)nH ^' (57) 

valid for both charge and flux qubits. Analytically we may proceed in the 
Breit-Wigner approximation of the transfer functions Xg/ei^)- We want to 
compare this dephasing rate with other dephasing mechanism affecting the 
qubit. For this purpose it is sufficient to make the approximation \x g \ ~ 
|Xe|; an d ignore effects of the quantum capacitance. Given that the rest 
of the integrand varies slowly compared to x(w) we arrive at the following 
expression for the dephasing rate for the charge qubit 

C Q 

= u Q Q—^-{n(u Q ) + l)n(u> ), (58) 

^osc 

where loq = 1/ \J LC osc and Q = 1/ '(ZqC^Lloq) is the approximate resonance 
frequency and Q value. A similar expression was found by Bertet et. al. 
considering a flux qubit coupled to a SQUID, in the high Q regime. 26 As 
expected the dephasing rate scales badly with the temperature due to the 
bosonic nature of the transmission line excitations. However we see that 
a low oscillator Q-value protects the qubit from dephasing, basically by 
pushing the noise spectrum up to higher frequencies. 



8. QUANTUM EFFICIENCY 

There is always a trade-off between speed in read-out and how fast the 
qubit is dephased. It is impossible to obtain information about the state of 
the qubit faster than the phase coherence between the states are lost. Given 
a longitudinal coupling between qubit and oscillator and a probing field in a 
coherent Glauber state, which is narrow in frequency distribution, we have 
derived the following general expressions for the dephasing and measurement 
rate 




where g is the coupling between qubit and oscillator degrees of freedom and 
x(w) describes the relation between charge (phase) of the oscillator and the 
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in field. We now assume that the function x(w) can be written on Breit- 
Wigner form as in Eq. (j35j) . Using this, we get explicit expressions for the 
dephasing rate in terms of circuit parameters 



( 



9\2 



(x e o? 



1 + 4 



1 + 4 



(60) 



and the phase difference between the two qubit states gets the simple form 

[2Q e 



2 arctan — a~[w ~ w n) 
V ^5 



2 arctan 



■(w-wg) 



(61) 



Choosing the drive frequency to be centered between the two resonance 
frequencies 



Y- 2 -, these expressions simplify to 

{gZo? 



r 



(xg) 2 



+ 



(xg) 2 



1 2 



^9^9 "^q) %ei%e + ^e) 

4^,, sin 2 ( arctan x e + arctan x 



v 9 I ' 



(62) 



where x g / e = Q 9 ^ £ (ujq — Uq)/u>q . Defining the quantum efficiency as rj = 
r ms /(2r^,) and using Eq. (|62|) the efficiency is given in terms of general pa- 
rameters in the Breit-Wigner approximation. Using the parameters given for 
the charge and phase qubit in Eqs. (|35|) and (|37j) respectively and neglect- 
ing second order effects in the qubit-oscillator coupling (small parameter is 
Cq/C% sc and Lq/L 2 respectively) we get a quantum efficiency rj = 1 for 
both charge and flux qubit. 

For another choice of drive frequency, we find an efficiency below unity. 
For these non-optimal frequencies, the charge magnitude on the oscillator 
differ for the two different qubit states. In this case, since the time to build 
up a charge on the oscillator depend on the qubit state, the time response 
of the oscillator will differ. Thus some information of the qubit state will be 
encoded in the response time, and since our setup only detects the phase, 
it misses this information, making the efficiency decrease. For the optimal 
drive frequency above, the charge on the oscillator is the same regardless of 
qubit state, and all information is encoded in the signal phase. 

In a more realistic situation, the qubit will of course be subject to var- 
ious additional sources of decoherence, i.e. thermal photons in the read-out 
lines and charge fluctuations in the environment, and these will reduce the 
quantum efficiency. Taking these sources into account is however beyond the 
scope of this paper. 
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9. MEASURING THE QUANTUM CAPACITANCE 
WITHOUT OSCILLATOR 

The calculation of the quantum efficiency for the charge and flux qubit 
was done assuming a rather sharp resonance in the function approxi- 
mating it with a Breit-Wigner form. We can however show that the result 
still holds if we remove the oscillator, thereby taking away the resonance 
completely. For this purpose we consider the circuit in Fig. |21 where the 
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Fig. 3. Cooper-pair box operated as a charge qubit. The read-out circuit 
only consists of a measurement capacitance C, still giving a full quantum 
efficiency rj = 1. 



charge qubit is directly coupled to the transmission line via the capacitor C. 
For this circuit the Hamiltonian is given by 

where C g / e = C/(l + gCa z ) and g = (8Eq) / \e 2 Ej) is the coupling between 
qubit and transmission line. The functions relating the in-field with the 
charge at the transmission-line boundary and the reflected signal are 

iZ c {uj)+Zo iZ c {uJ) + Z 

(64) 

where Zc{oj) = l/(wCg/ e ) is the impedance of the effective qubit-capacitor 
circuit. As in Sec. ISD we use a narrow bandwidth drive signal with drive 
frequency uo d an d Eq. (|59[) to calculate the dephasing and measurement rate 

r — r w ~ u/ I " r " 9> v 



(gZ ) 2 / A{uo d C g ) 2 4(oo d C e ) 2 
i,n 8u; 2 { l + (co d C g Z ) 2 + \ + {uj d C e Z Q ) 2 



J d 

t^ s = 4Ti n sin 2 ( arctan (uj d ZoC e ) — arctan (uj d ZoCg) ) . (65) 
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In this case a quantum efficiency of unity is also obtained, neglecting second 
order effects in the qubit-transmission line coupling g. This is however of 
more academic interest considering that typical parameters Zq = 50 Ohm, 
C g — C e = 1 fF, give a very small phase-difference on the order of Sip ~ 10~ 4 , 
even for a high drive frequency of aid = 1 GHz. A too small phase difference 
implies a small signal, which in a real setup will drown in the amplifier noise. 



10. CONCLUSION 

We have described dispersive charge and flux qubit readout using a slow 
oscillator as a quantum measurement process. Our main approximation is 
that the qubit follows the oscillator adiabatically, and we give conditions for 
this approximation to hold. Within this approximation the qubit acts as an 
effective linear circuit element in the oscillator circuit. For the charge qubit 
it is a capacitance for the flux qubit an inductance. The quantum dynamics 
of the effectively linear oscillator-transmission line is readily solved, and we 
get expression for all quantities in terms of the input field of the transmission 
line. Especially we treat the coupling between the oscillator and transmission 
line nonperturbatively. 

We use this formalism to calculate the quantum efficiency of the readout. 
This is found to be unity, within our approximation, for both charge and flux 
qubit readout. As an extra example we also find unit quantum efficiency for 
a charge qubit directly coupled to a transmission line. 



Appendix A. THE LEGENDRE TRANSFORM 

Given the Lagrangian as presented in Eq. and neglecting the po- 
tential part, 

L = -(t>iC ik (t>k - ak<fik, (A.l) 
we define the conjugated momenta according to 

dL 

Pj = —y = faCij - dj. (A.2) 

Given this, the Hamiltonian is obtained by a Legendre transformation of the 
Lagrangian 

H = pjfpj — L = pj(p k + a fe )(C" 1 )/ c j - 
\ 

- -(pi + Oj) (C~ )ijC jk (C~ 1 ) km (p m + a m ) + a.i{p k + afc)(C" 1 ) fe i = 

1 V v ' 

Si k 
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= PjiPk + afe)(C ,_1 )fc i - -{pi + a i )(C' 1 ) im (p m + a m ) + 
+ ai(p k + a fc )((7 _1 )fej = (jpi + a,i)(pk + a fe )(C _1 ) fe j - 

- ^fe + a *)(C -1 W(Pm + Om) = 2^ + a >)(C _1 )fa(Pm + «m), 

(A.3) 

making the derivation of the Hamiltonian a matter of inverting the matrix 
C, and subtracting the potential part of the Lagrangian. 



Appendix B. BREIT-WIGNER APPROXIMATION; 

CHARGE QUBIT + OSCILLATOR 

We are interested in the resonant behaviour of the transfer/response 
functions around u> = 1/ \] LC osc . We expand the relevant solution to N 9 ^ e (uj) = 
to first order in Zq 



UJ r 
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g/e 
23 
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leading to 
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, UJ — UJ, 
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2Q9/e J ~ "Qa/eC, 

Thus we approximate the transconductance as (assuming C c Zquj$ <C 1) 

-2 



g/e 



g/e 
Xo 



and the reflected phase as 
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